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ABSTRACT 

We study the generation of a stochastic gravitational wave (GW) background produced 
from a population of core-collapse supernovae, which form black holes in scenarios of 
structure formation. We obtain, for example, that the formation of a population (Pop- 
ulation III) of black holes, in cold dark matter scenarios, could generate a stochastic 
GW background with a maximum amplitude of /ibg — 10~ 24 and corresponding 
closure energy density of f^cw ~ 10~ 7 , in the frequency band fobs — 30 — 470 Hz 
(assuming a maximum efficiency of generation of GWs, namely, £GW max = 7 x 10 -4 ) 
for stars forming at redshifts z ~ 30 — 10. We show that it will be possible in the 
future to detect this isotropic GW background by correlating signals of a pair of 'ad- 
vanced' LIGO observatories (LIGO III) at a signal-to-noise ratio of ~ 40. We discuss 
what astrophysical information could be obtained from a positive (or even a negative) 
detection of such a GW background generated in scenarios such as those studied here. 
One of them is the possibility of obtaining the initial and final redshifts of the emission 
period from the observed spectrum of GWs. 

Key words: black hole physics — gravitation — cosmology: theory. 
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1 INTRODUCTION 

The detection of gravitational waves (GWs) will open up a 
new era in the history of astronomy and transform research 
in general relativity into an observational/theoretical study 
(Schutz 1999). 

The detection of GWs will directly verify the predictions 
of general relativity theory concerning the existence or not 
of such waves, as well as other theories of gravity (Thorne 
1987). The information provided by such waves is completely 
different when compared to that provided by electromag- 
netic waves. GWs carry detailed information on the coherent 
bulk motions of matter, such as in collapsing stellar cores or 
coherent vibrations of space-time curvature as produced, for 
example, by black holes. On the other hand, electromagnetic 
waves are usually an incoherent superposition of emissions 
from individual atoms, molecules and charged particles. 

There is a host of possible astrophysical sources of GWs: 
namely, supernovae, the collapse of a star or star cluster 
to form a black hole, inspiral and coalescence of compact 
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binaries, the fall of stars and black holes into supermassive 
black holes, rotating neutron stars, ordinary binary stars, 
relics of the big bang, vibrating or colliding of monopoles, 
cosmic strings and cosmic bubbles, among others (see, e.g., 
Thorne 1987, 1996, 1997; Schutz 1996, 1999). 

From the theoretical point of view there has been a 
great effort to study which are the most promising sources of 
GWs to be detected. In particular, the waveforms, the char- 
acteristic frequencies and the number of sources per year 
that one expects to observe are questions that have been 
addressed (see Thorne 1997; Schultz 1999; Grishchuk et al. 
2000 for a review). In a few years, starting from the observa- 
tions (waveforms, amplitudes, polarizations, etc.), it will be 
possible really to understand how GW emission is generated 
by the astrophysical sources. 

Because of the fact that GWs are produced by a large 
variety of astrophysical sources and cosmological phenom- 
ena, it is quite probable that the Universe is pervaded by 
a background of such waves. A variety of binary stars (or- 
dinary, compact or combinations of them), Population III 
stars, phase transitions in the early Universe and cosmic 
strings are examples of sources that could produce such pu- 
tative GW background (Thorne 1987). 

As GWs possess a very weak interaction with matter, 
passing through it without being disturbed, once detected 
they can provide information on the physical conditions from 
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the era in which they were produced. In principle, it will be 
possible to get information from the epoch when the galaxies 
and stars started to form and evolve. 

Concerning the production of GW backgrounds, it is 
worth mentioning that recently Blair & Ju (1996) and Fer- 
rari, Matarrese & Schneider (1999a,b) studied the cosmo- 
logical GW background produced by supernovae explosions 
that took place in the redshift range < z < 5. 

On the other hand, from considerations based on the 
Gunn- Peterson effect (Gunn & Peterson 1965), it is widely 
accepted that the Universe underwent a reheating (or reion- 
ization) phase between the standard recombination epoch 
(at z ~ 1000) and z > 5 (see Haiman & Loeb 1997; Loeb 
& Barkana 2001 for a review). However, at what redshift 
the reionization occurred is still an open question (Loeb 
& Barkana 2001), although recent studies conclude that it 
occurred at redshifts in the range 6 < z < 30 (Venkate- 
san 2000; Tegmark & Zaldarriaga 2000; Schmalzing et al. 
2000; Loeb & Barkana 2001). It is worth noting that present 
and future cosmic background radiation (CBR) studies can 
impose some constraints on the reionization phase of the 
Universe (see, e.g., Tegmark & Zaldarriaga 2000; Loeb & 
Barkana 2001). 

Although different models could possibly explain the 
reionization of the Universe, it is widely accepted that most 
of the contribution to the reionization is related to the 
formation and evolution of pre-galactic objects (sometimes 
called Population III objects) at high redshift (z > 10), such 
as subgalactic objects (M < 1O 9 M0) and stars formed from 
them (Loeb & Barkana 2001). This putative epoch where 
the formation of the Population III objects took place, and 
where the consequent reionization and reheating of the Uni- 
verse occurred, marked the end of an epoch named 'dark 
age' (see, e.g., Rees 1998). Also, note that the metalicity of 
~ W~ 2 Zq found in high — z Lya forest clouds (Songaila & 
Cowie 1996; Ellison et al. 2000) is consistent with a stellar 
population formed at z > 5 (Venkatesan 2000). 

The history of the Universe during the formation of the 
Population III objects could be investigated in the future 
with the New Generation of Space Telescope (NGST; Rees 
1998) and also, in principle, with GW observatories. Besides 
the reheating and reionization phases, putative Population 
III stars could produce GWs, particularly, from the forma- 
tion of neutron stars and black holes. Also, after the dark 
age, the epoch of the first light could be studied with large 
radio telescopes such as the Giant Meter Wavelength Tele- 
scope (GMRT) or the Square Kilometer Array Radio Tele- 
scope (SQA; Meiksin 1999). 

When the high-mass stars died as supernovae, they left 
stellar black holes as remnants. The formation of these stel- 
lar Population III black holes can, in principle, produce a 
GW background detectable by GW observatories. It is worth 
mentioning that a significant amount of GWs can also be 
produced during the formation of neutron stars. However, 
because this depends on the equation of state for the neu- 
tron star, which is not well defined, we consider here only the 
contribution of the black holes. Another possibility would 
be the generation of GWs through the so-called r-mode in- 
stability (Anderson 1998), which should be important for 
young, hot and rapidly rotating neutron stars, but we leave 
this issue for another study, to appear elsewhere. 

In the present study we have adopted a stellar gener- 



ation with a Salpeter initial mass function (IMF) as well 
as different stellar formation epochs. We then discuss what 
conclusions would be drawn whether (or not) the stochastic 
background studied here is detected by the forthcoming GW 
observatories such as LIGO and VIRGO. 

The paper is organized as follows. In Section 2 we 
present the basic ideas on the collapse of the first clouds and 
on the resulting stellar formation; in Section 3 we describe 
how to calculate the GW background produced during the 
formation of black holes in this scenario; in Section 4 wc 
present and discuss the numerical results; in Section 5 we 
consider the detectability of the putative GW background 
produced by the Population III black holes; and finally in 
Section 6 we present the conclusions. 



2 POPULATION III OBJECTS AND THE 
FIRST STARS 

The current theory of structure formation, based on cold 
dark matter (CDM) models, predicts that the first objects to 
collapse, the so-called Population III objects or mini-haloes, 
had a total mass of ~ 10 6 M Q and a formation epoch z ~ 
10—50 (see Tegmark et al. 1997, and references therein). The 
first stars, the Population III stars, started forming in these 
Population III objects of ~ 10 6 M Q and subsequently in more 
massive mini-haloes (see Haiman & Loeb 1997; Venkatesan 
2000; Loeb & Barkana 2001). 

As a result of the formation of the first stars, a pop- 
ulation of stellar black holes is formed after the supernova 
explosions associated with the high-mass stars. Then, know- 
ing the law of distribution of stellar masses, it is possible to 
obtain the number of stars that explode as supernovae, and 
so it is possible to determine the number (or event rate) of 
stellar black holes left as remnants. 

Thus, to proceed, the distribution function of stellar 
masses, the stellar IMF, for the first stars is required. Here 
the Salpeter IMF is adopted, namely 

<f>{m) = A m - (1+X \ (1) 

where A is the normalization constant and x = 1.35 (our 
fiducial value). The normalization of the IMF is obtained 
through the relation 

/•m a 

/ m4>(m)dm — 1, (2) 

J m\ 

where we consider mi = 0.1 M Q and m u = 125 M Q . It is 
worth noting that some authors argue (see, e.g., Gilmore 
2001) that there is evidence supporting the universality of 
the IMF, even for the first stars; on the other hand, other 
authors (see, e.g., Scalo 1998, among others) argue that the 
IMF may not be universal. In particular, the universality 
of the Salpeter exponent (x — 1.35) has been studied by 
recent evolutionary models for the Magellanic Clouds (de 
Freitas-Pacheco 1998). Some models, particularly those of 
the Large Magellanic Cloud, take into account constraints 
on the star formation history imposed by recent data on 
color-magnitude diagrams of field star clouds, showing that 
a steeper exponent, x = 2.0 is necessary to resolve the exces- 
sive production of iron obtained if one takes into account the 
Salpeter law (x = 1.35). Furthermore, concerning the star 
formation at high redshift, the IMF could be biased toward 
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high-mass stars, when compared to the solar neighborhood 
IMF, as a result of the absence of metals (Bromm, Coppi & 
Larson 1999, 2001). 

Then, for the standard IMF, the mass fraction of black 
holes produced as remnants of the stellar evolution is 



/bh 



M r <f>(m)dm, 



(3) 



where m m i n is the minimum stellar mass capable of produc- 
ing a black hole at the end of its life, and M r is the mass 
of the remnant black hole. Timmes, Woosley & Wheaver 
(1995) (see also Woosley & Timmes 1996) obtain, from 
stellar evolution calculations, that the minimal progenitor 
mass to form black holes is 18 < m m i„/Mg < 30 depend- 
ing on the stellar iron core mass. Thus, we assume that 
the minimum mass capable of forming a remnant black 
hole is m m i n = 25 Mq. For the remnant, M r , we take 
M r = am, where m is the mass of the progenitor star and 
q = 0.1 (see, e.g., Ferrari et al. 1999a, b). With these con- 
siderations at hand, the mass fraction of black holes reads 



/bh = 6.8 x 10- 



6.8 x 10" d for x = 1.35. 



To assess the role of possible IMF variations in our re- 
sults, other values of x have also been considered. Besides 
the standard IMF, two others have been studied, namely, 
with x — 0.3 and x = 1.85, which yield ten times and one- 
tenth of the mass fraction of black holes of the standard 
IMF, respectively. 

It is worth mentioning that stars formed with masses 
greater than 8 Mq to ~ 25 Mq also finish their lives as 
supernovae. Numerical studies have shown that these stars 
leave neutron star remnants, after forming iron cores with 
masses near the Chandrasekhar limit (Woosley & Timmes 
1996). These stars are important for injecting energy into 
the ambient medium and regulating the feedback of stellar 
formation. In the present paper only the generation of GWs 
that come from black holes formation have been studied, 
and so the progenitors of interest are stars with masses in 
the interval 25 < m/M < 125. 



3 GRAVITATIONAL WAVE PRODUCTION 

The GWs can be characterized by their dimensionless am- 
plitude, h, and frequency, v. The spectral energy density, the 
flux of GWs, received on Earth, F„, in erg cm~ 2 s _1 Hz _1 , 
is (see, e.g., Douglass & Braginsky 1979; Hils, Bender & 
Webbink 1990) 



16ttG ' 



(4) 



where u> bs = 27ri/ bs, with v bs the GW frequency (Hz) ob- 
served on Earth, c is the velocity of light, G is the gravita- 
tional constant and yfs^ is the strain amplitude of the GW 
(Hz" 1 / 2 ). 

The stochastic GW background produced by gravita- 
tional collapses that lead to black holes would have a spec- 
tral density of the flux of GWs and strain amplitude also 
related to the above equation (^). Therefore, in the above 
equation the strain amplitude takes into account the star 
formation history occurring at the 'first light', just after the 
'dark age' epoch. The strain amplitude at a given frequency, 



at the present time, is a contribution of black holes with dif- 
ferent masses at different redshifts. Thus, the ensemble of 
black holes formed produces a background whose character- 
istic strain amplitude at the present time is y^h- 

On the other hand, the spectral density of the flux can 
be written as (Ferrari et al. 1999a,b) 



F v = 



/^(fobs)d-RBH(m, z), 



(5) 



where f v {v bs) is the energy flux per unit of frequency (in 
erg cm- 2 Hz _1 ) produced by the formation of a unique black 
hole and dRsn is the differential rate of black holes forma- 
tion. 

The above equation takes into account the contribution 
of different masses that collapse to form black holes occur- 
ring between redshifts z c i and z c f (beginning and end of the 
star formation phase, respectively) that produce a signal at 
the same frequency v obs . On the other hand, we can write 
fviyobs) (Carr 1980) as 



(6) 



where /ibh is the dimensionless amplitude produced by the 
collapse to a black hole of a given star with mass m that 
generates at the present time a signal with frequency v b s - 
Then, the resulting equation for the spectral density of the 
flux is 



From the above equations we obtain for the strain amplitude 

(8) 



' = —J— I l'i:\\dli 



Thus, the dimensionless amplitude reads 



h%Q — / h%ndR, 



(9) 



(see de Araujo, Miranda & Aguiar 2000 for details). 

The dimensionless amplitude produced by the collapse 
of a star, or star cluster, to form a black hole is (Thorne 
1987) 



, 15 \ 1/2 GM r 

= 27T £GW ~ 



7.4 x io- 2 ° e y 2 l*^L w 



-gw 



M I \ IMpc 



(10) 



where £gw is the efficiency of generation of GWs and r z is 
the distance to the source. 

The collapse of a star to a black hole produces a signal 
with frequency (Thorne 1987) 



1/0,58 = 5^KG (1 + Z) 



i-3xio* Hz 0|)(i+*r\ 



(11) 



where the factor (1 + z) 1 takes into account the redshift 
effect on the emission frequency, that is, a signal emitted at 
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frequency v c at redshift z is observed at frequency i/ ODS 
u e (l + z)" 1 . The observed signal is in the range 



1.04 x i(r 

(1 + *=0 



Hz < Vobs < 



5.2 x 1(T 

(1 + z c! ) 



Hz, 



(12) 



obtained using the mass upper limit m u = 125 Mg , the mass 
lower limit m m i n = 25 Mq, and a = 0.1. 

For the differential rate of black hole formation we have 



gLR bh = p\(z)^—(f>(m)dmdz, 



(13) 



where p*(z) is the star formation rate (SFR) density (in 
Mq yr _1 Mpc -3 ) and dV is the comoving volume element. 

From the above equations we obtain for the dimension- 
less amplitude 



, a (7.4 x 10"' u a) 2 £GW 



fobs 



rn \ I flL 



M Q J V IMpc 



x ——<p(m)dmdz 

dz 



(14) 



In equation ( |l4| ) g(l is the luminosity distance to the 
source. The comoving volume element is given by 



dV = 4n y— ) r^^M, Ha, z)dz, 
with 



+ z) 2 (l + Hmz) - z{2 + z)Q A 
and the comoving distance, r z , is 

c „/ / ,„ , f z dz' 



where 

Qm = Odm + Ob and 1 = Ok + Am + Oa 



(15) 



(16) 



(17) 



(18) 



are the usual density parameters for the matter (M), i.e., 
dark matter (DM) plus baryonic matter (B), curvature (k) 
and cosmological constant (A). The function S is given by 

{sin a; if closed, 
x if flat, (19) 

sinhx if open. 

The comoving distance is related to the luminosity distance 
by 



<2l = r z (l + z). 



(20) 



The set of equations presented above can be used to find 
the dimensionless amplitude of the GW background gener- 
ated by black hole formation as a function of the SFR den- 
sity, and related to the 'first light' epoch. 

It is worth mentioning that the formulation used here 
is similar to that used by Ferrari et al. (1999a), but instead 
of using an average energy flux taken from Stark & Piran 
(1986), who simulated the axisymmetric collapse of a rotat- 
ing polytropic star to a black hole, we use equation ( |l(i| ) to 
obtain the energy flux, which takes into account the most 



relevant quasi-normal modes of a rotating black hole and 
represents a kind of average over the rotational parameter 
(see de Araujo et al. 2000). Both formulations present simi- 
lar results, since in the end the most important contributions 
to the energy flux come from the quasi-normal modes of the 
black holes formed, which account for most of the gravita- 
tional radiation produced during the collapse process. 

The SFR density, however, for the formation of the first 
stars is unknown. Star formation in other media could in 
principle give us some information on how things occurred 
in the 'first light' epoch, but unfortunately as one can see 
in what follows there are no compelling arguments in this 
direction. 

The formation of a bound cluster of stars requires a star 
formation efficiency of ~ 50 per cent when the cloud disrup- 
tion is sudden and ~ 20 per cent when cloud disruption 
takes place on a longer time-scale (Margulis & Lada 1983; 
Mathieu 1983; Ciardi et al. 2000). We define 'star formation 
efficiency' as the fraction of gas of a cloud that is converted 
into stars [this definition is similar to that used by Ciardi et 
al. (2000), among others]. 

On the other hand, Pandey, Paliwal & Mahra (1990) 
have investigated the influence of the IMF on the star for- 
mation efficiency, for clouds of different masses, and have 
concluded that the efficiency decreases if massive stars (the 
most destructive ones) are formed earlier. Some studies have 
also analysed the molecular gas properties and star forma- 
tion in nearby nuclear starburst galaxies (see, e.g., Planesas, 
Colina & Perez-Olea 1997) indicating the existence of giant 
molecular clouds with masses ~ 10 s - 1O 9 M , in which star 
formation process occurs in a short time (< 3 x 10 7 yr) with 
efficiency of conversion of gas into stars "Z 10 per cent. All 
these studies show the large uncertainties in the star forma- 
tion efficiency. 

Probably, the best way to infer the SFR density is to 
relate it to studies concerning the reionization of the Uni- 
verse. In the introduction of the present paper we argued 
that there are compelling arguments in favor of a reioniza- 
tion phase of the Universe, which probably occurred at red- 
shifts in the range 6 < z < 30 (Tegmark, Silk & Blanchard 
1994; Venkatesan 2000). There are at least two compelling 
reasons for reionization through Population III stars to be 
considered. First, the Population III stars are expected to 
form at z ^10, being capable of ionizing hydrogen. Secondly, 
the first stars create heavy elements, and can account for 
the metalicity of ~ 10~ 2 Zq found in Ly — a forest clouds 
(Venkatesan 2000). It is found that the amount of baryons 
necessary to participate in early star formation, to account 
for the reionization, would amount to a small fraction, /*, 
of all baryons of the Universe (see, e.g., Venkatesan 2000; 
Loeb & Barkana 2001). 

The above discussion suggests that we can write the 
SFR density as follows: 



dp* 
~~dt 



[n* Pc (i + z)% 



(21) 



where the term in brackets represents the stellar mass den- 
sity at redshift z, with p c the present critical density and 
f2* the stellar density parameter. The latter can be writ- 
ten as a fraction of the baryonic density parameter, namely, 
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Another relevant physical quantity associated with the 
GW background, produced by the first stars, is the clo- 
sure energy density per logarithmic frequency span, which 
is given by 

1 rfpcw 

"gw = r; ■ vi& ) 

p c a log u obs 

The above equation can be rewritten as 

^obs .-, 4-7T 2 2 ,2 /no\ 

"gw = — — r v = —jj2V obs h BG . (26) 

C p c ^IIq 

In the next section we present the numerical results and 
discussions, which come mainly from equation ( fl4| ) . Looking 
at this equation one notes that, to integrate it, one needs to 
choose the IMF, and to set values for the following param- 
eters: z ci , z c f, a, £gw, /*, H , f2 B , fi M and Qa- 



4 NUMERICAL RESULTS AND DISCUSSIONS 

Based on models of structure formation, the first objects 
that collapsed should have had masses around 10 6 Mq (see, 
e.g., Tegmark et al. 1997, among others), and objects with 
higher masses should have collapsed subsequently. Evi- 
dently, as the density fluctuations could have had peaks 
higher than la values, clouds could have collapsed earlier. 
As a result, there were clouds with different masses collaps- 
ing around the redshift of collapse of 10 6 Mq. 

We have considered that the first stellar formation is re- 
lated to the collapse of Population III objects. To evaluate 
the GW background produced by the formation of the Pop- 
ulation III black holes, it is necessary to know the redshifts 
at which they began and finished being formed. This is a 
very hard question to answer, since it involves knowledge of 
the role of the negative and positive feedbacks of star forma- 
tion, which are regulated by cooling and injection of energy 
processes. 

Should the stochastic GW background studied here be 
significantly produced and detected at a reasonable confi- 
dence level, the present study can be used to obtain the red- 
shift range where the Population III black holes were formed, 
independently of any CDM modelling. In Fig. 1 an exam- 
ple is given of how one could get z c i and z c t from the curve 
/ibg versus v obs . Knowing the frequency band 

^min £Anax 

and using equation (|l2j), one obtains z c \ and z c { (see Fig. 1), 
which are, therefore observable. Note that we have assumed 
as did Ferrari et al. (1999a, b) that a is a constant (a = 0.1). 

Is is worth noting that a may depend sensitively on 
the metallicity: the lower the value of Z, the higher are the 
remnant masses and the less ejected material there is relative 
to Zq stars. More realistically there would be a dependence 
of a on the progenitor mass. On the other hand, the value 
a — 0.1 adopted can be considered as a mean value for 
the progenitor masses studied here. If a is not very well 
determined, this would mean that the observed frequency 
band does not uniquely fix the redshift band where the black 
holes are formed. 

In the introduction of the present paper we mention 
that different studies related to the reionization of the 
Universe set this epoch as being somewhere in the range 
6 < z < 30. It is also clear that the reionization process is 
not instantaneous. Stars start forming at different redshifts, 



Z 




10 100 1000 



obs 

Figure 1. Example of how one could obtain from /ibg versus f b s 
the initial (final) redshift z c \ (z c f) of the GW emission period. We 
have adopted a = 0.1. 



creating ionized bubbles (Stromgren spheres) around them- 
selves, which expand into the intergalactic medium (IGM), 
at a rate dictated by the source luminosity and the back- 
ground IGM density (Loeb & Barkana 2001). The reioniza- 
tion is complete when the bubbles overlap to fill the entire 
Universe. Thus the epoch of reionization is not the epoch of 
star formation. There is a non-negligible time-span between 
them. Here, we have chosen different formation epochs to see 
their influence on the putative GW background and also to 
see if it could be detected by the forthcoming GW antennas. 

The first relevant quantity appearing in the equation for 
the GW background is the IMF. As discussed in the previous 
section, we have adopted the standard IMF [x = 1.35), as 
our fiducial case, and have also studied two other cases to 
assess the role of the IMF variations on /ibg- 

The second relevant parameter is egw, the efficiency 
of production of GWs, whose distribution function is un- 
known. Thus, we have parameterized our results in terms of 
its maximum value, namely, £Gw max = 7 x 10~ 4 , which is 
obtained from studies by Stark & Piran (1986) who simu- 
lated the axisymmetric collapse of a rotating star to a black 
hole. We will see below that, if £gw is a very tiny fraction of 
the maximum value, the detection of the GW background 
whose existence we propose, is very improbable, even for 
advanced antennas. 
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To calculate /ibg we still need to know f2* , which has a 
key role in the definition of the SFR density. As discussed in 
the previous section studies related to the reionization of the 
Universe can shed some light on f2*. From different studies 
one can conclude that a few percent, maybe up to ~ 10 per 
cent, of the baryons must be condensed into stars in order for 
the reionization of the Universe to take place. Here we have 
set the value of f2* in such a way that it amounts to 1 per 
cent of all baryons (our fiducial value). In the next section 
we discuss the detectability of the GW background, and we 
then parameterize our results in terms of /* = fi^/fie- 

Looking at equation ( |l4| ) one could think it would de- 
pend critically on the cosmological parameters: Ho, Ob, 
J7dm, Oa- However, our results show that, given the red- 
shifts involved in our calculations, /ibg depends only on Ho 
and Qb', the latter dependence occurs because this param- 
eter appears in the SFR density. The quantity hi 00 fls = 
0.019 ± 0.0024 (where /iioo is the Hubble parameter given 
in terms of 100 kms -1 Mpc -1 ) is obtained from big bang 
nucleosynthesis studies (see, e.g., Buries et al. 1999). 

In Table 1 we present the redshift band, z c i and z c f for 
the models studied and the corresponding GW frequency 
bands. For the cosmological parameters we have adopted 
h wo = 0.65, n M = 0.3, fl B = 0.045 and fi A = 0.7. Keep 
in mind that our results are sensitive to the combination 
/iioo^B- We have also adopted a — 0.1, /* = 0.01 and the 
standard IMF. 

Note that no structure formation model has been used 
to find the black hole formation epoch. Instead we have sim- 
ply chosen the values of z to see whether it is possible to 
obtain detectable GW signals. In the next section it will be 
seen that, unless £gw is negligible, the GW background that 
we propose here can be detected. Our choices, however, can 
be understood as follows. The greater the redshift forma- 
tion, the more power the masses related to the Population 
III objects have. Thus, of our models A to D, our model 
D (A) has more (less) power when compared to the oth- 
ers. Models E, F and G would mean a more extended star 
formation epoch, which means that the feedback processes 
of star formation are such that they allow a more extended 
star formation epoch when compared to models B, C and 
D, respectively. 

Concerning the reionization epoch, as already men- 
tioned, it occurred at lower redshifts as compared to the first 
star formation redshifts. Loeb & Barkana (2001) found, for 
example, that if the stars were formed at z ~ 10 — 30, with 
standard IMF, they could have reionized the Universe at 
redshift z ~ 6. Our models A, B and E, for example, could 
account for such a reionization redshift. 

In addition, we also consider a model with z c t = 5 (see 
model H of Table 1) to verify if a final epoch of star forma- 
tion close to this redshift could produce a detectable signal 
for the VIRGO and LIGO experiments. 

If the process of structure formation of the Universe and 
the consequent star formation were well known, one could 
obtain the redshift formation epoch of the first stars. On the 
other hand, as discussed below, if the GW background really 
exists and is detected, one can obtain information about the 
formation epoch of the first stars. 

A relevant question is whether the background we study 
here is continuous or not. The duty cycle indicates if the 
collective effect of the bursts of GWs generated during the 



Table 1. The redshifts of collapse for our models and the cor- 
responding GW frequency bands. The cosmological parameter 
hl 0Q Q B = 0.019 (see the text), a = 0.1, /* = 0.01 (our fidu- 
cial value) and the standard IMF are adopted. 



Model z c i z c { Au 
(Hz) 



A 


20 


10 


50-470 


13 


30 


20 


34-250 


C 


10 


30 


25-170 


D 


50 


40 


20-130 


E 


30 


10 


34-470 


F 


40 


10 


25-470 


G 


50 


10 


20-470 


H 


15 


5 


65-870 



collapse of a progenitor star generates a continuous back- 
ground. The duty cycle is defined as follows: 

DC = I oLRbhAtgw(1 + z), (24) 

2 cf 

where Atgw is the average time duration of single bursts 
at the emission, which is inversely proportional to the fre- 
quency of the lowest quasi-normal mode of the rotating 
black holes (see, e.g., Ferrari et al. 1999a), which amounts 
to ~ 1 ms for the mass range of the black holes considered 
here. 

Since the star formation rate could be high, a signifi- 
cant amount of GWs could be produced. We also note that, 
independently of the primordial cloud mass and of redshift 
of collapse, star formation occurring at high redshift could 
produce high duty cycle values, which lead us to conclude 
that the stochastic GW background could be continuous. 
For all the models studied here the duty cycle is 3> 1. 

The amplitude /ibg of the GW background in terms of 
egw (the efficiency of generation of GWs) , in the frequency 
band p-min — fmax is shown in Fig. 2 for the models A, E and 
F of Table 1. In the next section we discuss the detectability 
of such a putative background. 

Note that, the earlier that star formation occurs, the 
greater is the GW amplitude /ibg- This can be explained as 
follows. For a given individual source, the higher the red- 
shift, the lower is the amplitude of the GWs generated. On 
the other hand, the higher the redshift of star formation, 
the greater is the SFR density. As a result there is a more 
significant overlapping of bursts of GWs at higher redshifts. 

We find, for example, that the formation of Population 
III of black holes, in model D, could generate a stochastic 
GW background with amplitude /ibg ^ (0.8 - 2) x 10~ 24 
and a corresponding closure density of Ogw — (0.7 — 1.4) x 
10~ 8 , in the frequency band v a ^ s ~ 20 — 130 Hz (assuming 
an efficiency of generation egw — 7 x 10 -4 , the maximum 
one). 

It could be argued that the formation of stars could not 
be restricted to the epoch of the collapse of the first Popu- 
lation III objects. A new surge of star formation associated 
with the collapse of Population III objects of M > 1O 6 M0 
at lower redshifts could occur. The existence of this new 
surge of star formation would depend on the role of the neg- 
ative and positive feedbacks, which are regulated by cooling 
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Figure 2. The background amplitude of the GWs as a function 
of j/ b s and the efficiency of GW generation £gw for models A, 
E and F of Table 1 



and injection of energy processes of the previous star for- 
mation surge. If another surge of star formation took place, 
then another GW background could be generated, and a 
partial superposition with the background previously gener- 
ated could also occur. As a result for some frequency bands 
the GW amplitude could be enhanced. 

Another possibility would be a star formation process 
taking place during the time of collapse of the first Popula- 
tion III objects and continuing during the collapse of objects 
of higher masses. In such a case the star formation would 
occur for a large redshift span, and as a result the frequency 
band, the amplitude and the closure energy density of GWs 
could be larger. As before, the role of the negative and posi- 
tive feedbacks of the star formation would be the key point. 
In the models E, F and G we have considered such a possi- 
bility (see also Fig. 2). 

Certainly, the GW background produced depends on 
the star formation history. A different star formation his- 
tory would produce different results for both the values of 
/ibg and frequency bands. However, it is hard to avoid the 
conclusion that, if the first stars had been formed at high 
redshift, a significant amount of GWs would have been pro- 
duced as well. 

To assess the role of possible IMF variations we have 
considered, besides the standard IMF, two others, namely, 



with x = 0.3 and x = 1.85, which yield ten times and one- 
tenth of the mass fraction of black holes of the standard 
IMF, respectively. For the model A, our calculations show 
that for x = 0.3 (x = 1.85) the maximum /ibg is a factor 
~ 3 (~ 4) greater (smaller) than case with standard IMF 
(x = 1.35). In the next section we consider the role of the 
IMF variations on detectability of the background of the 
GW background that we propose exists. 



5 DETECTABILITY OF THE BACKGROUND 
OF GRAVITATIONAL WAVES 

The background predicted in the present study cannot be de- 
tected by single forthcoming interferometric detectors, such 
as VIRGO and LIGO (even by advanced ones). However, it 
is possible to correlate the signal of two or more detectors 
to detect the background that we propose exists. Michelson 
(1987) was the first to show that this kind of signal can, in 
principle, be detected by correlating the outputs of two dif- 
ferent detectors. However, the main requirement that must 
be fulfilled is that they must have independent noise. This 
study was improved by Christensen (1992) and by Flanagan 
(1993). The reader should also refer to the papers by Allen 
(1997) and Allen & Romano (1999) who also deal in detail 
with such an issue. 

To assess the detectability of a GW signal, one must 
evaluate the signal-to-noise ratio (S/N), which for a pair of 
interferometers is given by (see, e.g., Flanagan 1993; Allen 
1997) 



(S/N) 2 



( 9 A)t 



<<V)4 2) 



(25) 



where is the spectral noise density, T is the integration 
time and 7(1/) is the overlap reduction function, which de- 
pends on the relative positions and orientations of the two 
interferometers. For the 7(f) function we refer the reader to 
Flanagan (1993), who was the first to calculate a closed form 
for the LIGO observatories. Flanagan (1993; see also Allen 
1997) showed that the best window for detecting a signal is 
< v < 64 Hz, where the overlap reduction function has 
the greatest magnitude. 

Here we consider, in particular, the LIGO interferom- 
eters. Their spectral noise densities have been taken from 
a paper by Owen et al. (1998)- who in turn obtained them 
from Thorne, by means of private communication. 

In Table 2 we present the S/N for one year of ob- 
servation with a — 0.1, Qb/iioo = 0.019, /* = 0.01 and 
£Gw max = 7 x 10~ 4 for the models of Table 1, for the three 
LIGO interferometer configurations. 

Note that for the 'initial' LIGO (LIGO I) there is no 
hope of detecting the GW background we propose here, even 
for ideal orientation and locations of the interferometers, 
i.e., \y{y)\ = 1. For the 'enhanced' LIGO (LIGO II) there is 
some possibility of detecting the background, since S/N > 1, 
if £gw is around the maximum value. Even if the LIGO II 
interferometers cannot detect such a background, it will be 
possible to constrain the efficiency of GW production. 

The prospect for the detection with the 'advanced' 
LIGO (LIGO III) interferometers is much more optimistic, 
since the S/N for almost all models is significantly greater 
than unity. Only if the value of £gw were significantly lower 
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Table 2. For the models of Table 1 we present the S/N for pairs of 
LIGO I, II and III ('first', 'enhanced' and 'advanced', respectively) 
observatories for one year of observation. Note that an efficiency 
of generation £GW max = 7 X 10~ 4 is assumed. 



Model 


LIGO I 


S/N 
LIGO II 


LIGO III 


A 


8.3 x 10" 3 


1.6 


6.6 


B 


8.5 x icr 3 


2.3 


26 


C 


8.7 x icr 3 


2.7 


47 


D 


8.1 x icr 3 


2.5 


51 


E 


2.7 x icr 3 


5.7 


37 


F 


5.0 x icr 3 


12 


120 


G 


7.7 x 10- 2 


21 


260 


H 


4.6 x icr 3 


0.5 


1.7 



than the maximum value would the detection not be possi- 
ble. In fact, the S/N is critically dependent on this parame- 
ter, whose distribution function is unknown. 

Note, for example, that it is possible to detect a GW 
background with the 'advanced' LIGO, even for star for- 
mation for which z c f ~ 5 (model H) , if £gw is around the 
maximum value. 

Let us now look at how the variations of the parameters 
modify our results. First of all, note that the larger the star 
formation redshift band, the greater is the S/N. Secondly, 
the earlier the star formation, the greater is the S/N. It 
is worth recalling that, if one can obtain the curve of /ibg 
versus v ohs and if the value of a is known, one can find the 
redshift of star formation. 

The S/N is also sensitive to variations of a. The larger 
a, the lower are the GW frequencies and the higher is /ibg, 
and since the best window for detection is around < v < 
64Hz, the S/N is higher. 

Even if a is not known beforehand, it is possible to im- 
pose a constraint on its values, and also on the redshift star 
formation epoch. For example, if one found from GW ob- 
servations that the GW frequency band were 40 — 200 Hz, 
one would obtain (using Equation 11) that a ~ 0.1 — 0.4 and 
Z{ ~ 5 — - 50. On the other hand, if one knew that the star for- 
mation redshift band were Zf ~ 10 — 30, through some model 
of structure formation or whatever observational data, and 
the GW frequency band were known, say 40 — 200 Hz, one 
would obtain that a ~ 0.1. 

It would be interesting to perform a study considering 
q as a function of the progenitor mass, which would result 
in a more realistic model. There are some studies in the 
literature considering how the remnant mass depends on 
the progenitor (see, e.g., Fryer & Kalogera 2001), but we 
will not consider such an issue here. 

To assess the role of the variations in the IMF, we mod- 
ify its exponent x. For x = 0.3 (x — 1.85), the S/N is ~ 10 
(~ 0.1) times the S/N of the standard IMF. As expected for 
an IMF biased toward high- (low-) mass stars, where one has 
a greater (lower) amount of black holes, the S/N is greater 
(lower) . 

Note that the S/N for a given formation epoch, IMF 
and a, and for one year of observation, still presents a de- 
pendence on QbIi 2 , /* and £gw, namely 

S/N oc n B h 2 f* £gw- (26) 



The value of Qsh 2 is well constrained by primordial nucle- 
osynthesis studies. For /* we have adopted a value of 0.01, 
which is a very conservative choice. Note that the value 
for this parameter can be obtained from studies concern- 
ing the reionization of the Universe, which is very difficult 
to model, but there is some agreement in the literature (see, 
e.g, Gnedin 2000; Venkatesan 2000; Loeb & Barkana 2001) 
among the different models of the reionization of the Uni- 
verse which lead us to conclude that /* could range from a 
few up to 15 per cent. 

For £gw, the situation is more complicated since its 
distribution function is unknown. We have adopted here the 
maximum value as a reference, but if its actual value is much 
less than this value the S/N could be lower than unity for 
all the models studied here, even for a LIGO III pair. Let us 
think of what occurs with other compact objects, namely, 
the neutron stars, to see if we can learn something from 
them. Hot and rapidly rotating neutron stars can lose an- 
gular momentum to gravitational radiation via the so-called 
r-mode instability (Anderson 1998). This could explain why 
all known young neutron stars are relatively slow rotators. 
The black holes could have had a similar history, i.e., they 
could have been formed rapidly rotating and lost angular 
momentum to gravitation radiation via their quasi-normal 
modes. If this was the case, the value of £gw could be near 
the maximum one, or in the worst case it could have a value 
to produce S/N > 1 at least for a LIGO III pair. 

In order to assess the values of /* and £gw that yield 
S/N > 1, for a given formation epoch, IMF, a and fish 2 , we 
present in Fig. 3 the regions in the (/* , £gw) plane where 
S/N could be greater than unity for a pair of LIGO III in- 
terferometers. Note that unless £gw is very small, S/N can 
be significantly greater than unity, indicating that the back- 
ground could in principle be detected in the near future. 

A relevant issue is whether there are other GW back- 
grounds that could be confused with that of the present 
study. Relic GWs generated in the very early Universe can 
in principle present a signal in the LIGO bandwidth. The or- 
dinary inflationary models, however, predict Ogw ~ 10~ 15 
(see, e.g., Schultz 1999; Giovannini 2000; Maggiore 2000a,b). 
This is much less than our models predict, and therefore 
undetectable even with a pair of LIGO III interferometers. 
Other models, such as string cosmologies, provide different 
predictions (see, e.g., Schutz 1999) with values of Ogw that 
could be much greater than our studies predict, which could 
render the background studied here undetectable. Other 
GW backgrounds exist in the bandwidth of LIGO, which 
could have been produced at < z < 5, namely: (a) a cosmo- 
logical population of core-collapse supernovae (Ferrari et al. 
1999a); (b) a population of young rapidly rotating neutron 
stars (Owen et al. 1998, Ferrari et al. 1999b); and (c) dou- 
ble neutron star binaries (Schneider et al. 2000). The first 
two backgrounds, however, have energy density shifted for 
higher GW frequencies when compared to our predictions. 
Moreover, the S/Ns of these backgrounds are less than unity, 
even for a pair of LIGO III. The last one has a frequency 
band ranging from ~ 10~ 5 Hz up to ~ 10 2 Hz, and therefore 
there is a partial overlap with the background of our study. 
The GW amplitudes, however, would be comparable only if 
£Gw/£GW max << 1. 



GW background from black hole formation 9 



O 0.01 



GO 



0.001 




0.02 0.04 0.06 0.08 0.10 



Figure 3. Relative efficiency of GW generation, £Gw/£GW ma x> 
as a function of the fraction of baryons participating in the early 
star formation, /*, for models A, E and F of Table 1 for a pair of 
LIGO III interferometers. The curves represent where the S/N = 
1; above them S/N > 1. As in Table 2, E G w m „ = 7 X 10 -4 is 
assumed. 



6 CONCLUSIONS 

We present here a study concerning the generation of GWs 
produced from a cosmological population of black holes. 
These objects are formed as a consequence of the collapse of 
pre-galactic objects (Population III objects) that form the 
first generation of stars at high redshift. Our results show 
that different structure formation models which predict the 
formation of the first objects at z > 10 could, in principle, 
predict the formation of pre-galactic black holes and a sig- 
nificant stochastic GW background associated with them. 
Our results lead us to conclude that star formation occur- 
ring at high redshifts could have large duty cycles and so 
the stochastic GW background generated is continuous. 

We consider that stars are formed following a Salpeter 
IMF and having masses in the range 0.1 — 125M0 . Certainly, 
the res mlta prcacntcd here arc dependent on thio particular 
choice. l A steeper IMF would modify the number of high — 



mass stars, modifying the peak of /ibg and the frequency 
band of the GWs. For an IMF with x = 0.30 (x = 1.85) the 
IMF is biased toward high- (low-) mass stars, as a result the 
S/N is ~10 (~0.1) times the S/N predicted with the use of 
a standard IMF. It would be of interest, however, to have 



a look in detail at studies of the metallicity of high-z Ly a 
clouds to see if it is possible to constrain the Population III 
IMF. 

If we consider £gw — 7.0 x 10 -4 (see, e.g., Stark 
& Piran 1986) then we obtain h B a ^ (0.8 - 2) x 10" 24 and 
tt GW ~ (0.7-1.4) xlO -8 at u ohs ~ 20- 130 Hz for the model 
D. Thus, this GW background produced as a consequence of 
the formation of the first stars in the Universe is capable of 
being detected by a pair of 'advanced' LIGO interferometers. 

As seen, with reasonable parameters, our results show 
that a significant amount of GWs is produced related to the 
the Population III black hole formation at high redshift, and 
can in principle be detected by a pair of LIGO II (or most 
probably by a pair of LIGO III) interferometers. However, 
a relevant question should be considered: What astrophysi- 
cal information can one obtain from whether or not such a 
putative background is detected? 

First, let us consider a non-detection of the GW back- 
ground. The critical parameter to be constrained here is 
egw- A non-detection would mean that the efficiency of 
GWs during the formation of black holes is not high enough. 
Another possibility is that the first generation of stars is 
such that the black holes formed had masses > IOOMq , and 
should they form at z > 10 the GW frequency band would 
be out of the LIGO frequency band. 

Secondly, a detection of the background with a signif- 
icant S/N would permit us to obtain the curve /ibg versus 
fobs- From it, as discussed above, one can constrain a and 
the redshift formation epoch; and for a given IMF and Q, B h 2 , 
one can also constrain the values of /* and egw- On the 
other hand, using the curve /ibg versus z/ b s and in addition 
other astrophysical data, say CBR data, models of structure 
formation and reionization of Universe, a constraint on egw 
can also be imposed. 
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